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Abstract 

This paper concerns the Routh reduction procedure for Lagrangians systems with symmetry. 
-Y«, i It differs from the existing results on geometric Routh reduction in the fact that no regularity 

conditions on either the Lagrangian L or the momentum map Jl are required apart from the 
momentum being a regular value of Jl- The main results of this paper are: the description of 
a general Routh reduction procedure that preserves the Euler-Lagrange nature of the original 
system and the presentation of a presymplectic framework for Routh reduced systems. In 
addition, we provide a detailed description and interpretation of the Euler-Lagrange equa- 
tions for the reduced system. The proposed procedure includes Lagrangian systems with a 
non-positively definite kinetic energy metric. 
MSC2000: Primary 70H33. Secondary 70G65, 70H03, 53D20. 
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1 Introduction 

In geometric accounts to Routh reduction [U [HI [H] it is custom to first consider the restriction of 
the dynamics associated to a Lagrangian system with symmetry to the level set of the momentum 
map J£ (/i) and then reduce to the quotient manifold under the action of the isotropy group 
G M . Because this procedure is very similar to symplectic or Marsden-Weinstein reduction, one 
often states that Routh reduction is the Lagrangian analogue of cotangent bundle reduction (see 
e.g. [9J [H]). Another essential ingredient in the work on Routh reduction is the existence of a 
regularity condition on the momentum map Jl- it guarantees that the manifold J^ 1 (^)/G M is 
diffeomorphic to T(Q/G) x Q/G^. Roughly said, this condition provides the manifold J i ^ 1 (/z)/G /i 
with a tangent bundle structure, an essential feature to reinterpret the reduced dynamics as being 
Euler-Lagrange equations. We say that Lagrangians satisfying this regularity condition are G- 
regular. For example, the classical T — V Lagrangians fit in this context. Unfortunately, there 
are many Lagrangians where this procedure can not be carried out. This is the case for singular 
Lagrangians, though there are also instances of regular Lagrangians where it is impossible to 
perform Routhian reduction in this sense. A simple one is the Lagrangian L(q 1 ,q 2 ,q 1 , q 2 ) = 
(q 1 ) 2 + q^q 2 — V(q 1 ) defined in Q = K 2 and symmetric with respect to the action of the group 
G = K of translations along the q 2 coordinate (using the old language in Mechanics, q 2 is a cyclic 
coordinate). The momentum map is Jl = q . For a frozen value ^ — q 1 of Jl, the quotient 



(TQ,L) - (TQ,Rn 

MW 

{Jl\^/G^X) » (TQ/G^TZn = (t(Q/G) x (Q/G^ x ~g),TzA 

Figure 1: Schematic drawing of reduction on a tangent bundle 

J£ (/z)/G M can not be identified with T(Q/G) x Q/G^ = T(Q/G) as the coordinates of the 
quotient space J£ 1 {fj,)/G fl are precisely {q x ,q 2 ) whereas the coordinates oiT(Q/G) are {q 1 ,^ 1 ). 
See M3.2I below. The main result of this paper is the generalization of Routh reduction to a 
reduction technique that holds for arbitrary Lagrangians with the only requirement of [i being 
a regular value of Jl- The price paid for that is that the reduction process will be carried out 
in the entire space TQ and not only on J£ (/i). The fact that the reduction takes place in the 
whole space does not mean that we loose control of the momentum, at the contrary, the reduction 
process will keep track of the value of Jl in the same spirit of the classical Routh reduction. 

In order to relate the presented reduction technique to existing results on geometric Routh reduc- 
tion, we first mention that the correspondence between Routh reduction and Marsden-Weinstein 
reduction holds at a much more fundamental level: Routh-reduction is Marsden-Weinstein re- 
duction applied to the Poincare-Cartan symplectic structure on the tangent bundle (see [12]). 
This correspondence was then used to generalize Routh reduction to Lagrangians that are quasi- 
invariant (up to a total time derivative). 

Here we wish to generalize Routh reduction from a different point of view. In classical treatments 
on Routh reduction of a Lagrangian system with cyclic coordinates, the reduced system is again 
a Lagrangian system with a new Lagrangian which is called the Routhian. It is this observation 
what we take as a starting point in this paper: Routh reduction is a reduction technique that 
preserves the Euler-Lagrange nature of the original system. 

The different steps in the proposed reduction are best illustrated by means of the schematic 
diagram in Fig. [T] For that purpose, we shall denote for now a Lagrangian system as a couple 
{TQ,L), Q being the configuration manifold and L a function on TQ called the Lagrangian. 
The first step is to consider an equivalent Lagrangian system on Q, with Lagrangian R^. This 
new Lagrangian has the property that solutions of the Euler-Lagrange equations of {TQ, L) with 
momentum /i are solutions to {TQ, R^) with momentum 0. For the sake of completeness we 
mention here that {TQ, R 11 ) is not conservative and additional force terms should be taken into 
account; and that the function R^ is only G M invariant but, as will become clear, we can make 
sense to a momentum map of i? M taking values in q*, the dual of the Lie-algebra of G. This step 
is not new and was carried out, for instance, in |S]. The Routh reduction technique described 
in [9] is schematically presented in Fig. [1] by the arrows 1 and 2, followed in this order. 

We shall follow a different path: in our approach the second step is to perform a reduction on 
the system {TQ, _R M ) that is similar to Lagrange-Poincare reduction, i.e. we reduce the system 
to {TQ/G^, TZ^), where TV 1 is the quotient of i? M (represented by arrow 3 in the above diagram). 
We show that this step, in contrast to general Lagrange-Poincare reduction, preserves the Euler- 
Lagrange nature of the original system. To explain this statement in more detail, we mention 
here that TQ/G^ can be identified with the fibred product of three bundles over Q/G: T{Q/G), 
Q/Gfj, and § (the latter being the bundle associated to g). Roughly said, the Routhian W 1 



depends on velocities v x in T(Q/G) and on points y and £ in the fibres of Q/G^ and g above x 
respectively. When we say that lZ^(v x ,y,t;) is a Lagrangian, we understand that the variables y 
and £ are to be interpreted as configuration coordinates. This might be surprising, since £ is the 
projection of the vertical part of a velocity upstairs, i.e. a tangent vector in TQ. 

It is immediately clear now that the Lagrangian W* is singular, i.e. it does not depend on 
the velocities of the variables y and £. This is another price we pay for keeping the Euler- 
Lagrange nature of the system: we have to work with singular Lagrangians. But the singularity 
of these systems is of a specific type, which we call throughout this paper intrinsically constrained 
Lagrangian systems. They are studied in detail in Section [2j It is well known that singular 
Lagrangians contain constraints on the solutions to the associated Euler-Lagrange equations. It 
is shown that one of these intrinsic constraints embedded in the system (TQ/G^, TV 1 ) is a reduced 
version of the fixed momentum condition Jj, = \x. As a result of this Lagrangian interpretation of 
the reduced system (TQ/G^TV 1 ), we shall define a presymplectic formulation and, in the most 
general case, one can apply the presymplectic constraint algorithm [6l [7j [8] to find solutions to 
the Euler-Lagrange equations. 

Depending on the nature of these constraints we describe in more detail two distinct cases where 
this reduced singular Lagrangian system is equivalent to a variational problem on the space 
J z ^ 1 ( y u)/G /J ; the first case leads to standard Routh reduction and the second case is new and, 
for instance, is applicable to geodesies in general relativity where the metric is invariant under a 
lightlike vector field. 

Finally, we wish to mention that besides the importance of Routh reduction for Mechanical 
systems by itself, the special features of this process may shed light to other situations. In 
particular, this might be the case for Field Theories. The covariant reduction of these Lagrangian 
systems under the action of a group of symmetries has been developed in recent years (see for 
example |14])„ though the Hamiltonian counterpart is much less explored. A generalization of 
Routh techniques to this setting could be of much interest to tackle some of the difficulties 
encountered in this context (for instance, the lack of meaning of the notion of fixing the value of 
the momentum map). 

Throughout this paper manifolds are always assumed to be smooth finite dimensional (Hausdorff, 
second countable, C°°) and smooth always means of class C°° . We will often consider the pull- 
back of a bundle, a function, a one- or two-form along a map / : M — > N between two manifolds. 
To reduce the notational complexity in this paper, we sometimes denote the pull-backed object 
with the same symbol as the object itself. It should be clear from the context what is meant. 
Given two bundles ~k\ : M\ — > N and TT2 : Mq, — > N over a manifold N, then we often consider the 
fibred product bundle Mi x jy M.% , or simply M\ x M-i when no confusion is possible. Elements in 
such a fibred product are denoted as a couple (mi, mz), with mi £ M\ and 1112 £ -M2 such that 
TTi(mi) = ^2(7712). Finally, velocities at a point m £ M are typically denoted by v m , w m £ T m M. 
When we consider the tangent vector to a curve m : I — > M at a time t, we shall write m(t) for 
the corresponding curve in TM. 



2 Lagrangian systems on fibred manifolds 

We start with recalling some standard concepts in Lagrangian mechanics. Assume that L is a 
Lagrangian function defined on the tangent space of a manifold M . 

Definition 1 A Lagrangian system is a triple (M,L,F) with L a function on TM and F a 
co-vector valued function TM , i.e. F : TM — ► T*M, fibred over the identity. The map F is 



called the force term of the Lagrangian system. A Lagrangian system (M, L, F) is said to be 
conservative if F = 0. 

In the above notion of conservative systems we assume that conservative forces are always taken 
into account by means of their potential energy in the Lagrangian function L. 

Definition 2 A curve m : I = [a, b] — > M is said to be a critical curve for the Lagrangian system 
(M,L,F) if for arbitrary variations 5m with fixed endpoints 5m(a,b) = 0, 

5 J L(m(t))dt + I (F(m(t)),6m(t))dt = 

holds. 

It is well-known that a critical curve satisfies 

£C(L)(m(t))+F(m(t)) = 0, 

with £C the Euler-Lagrange operator £C{L) : T^'M — > T*M. When expressed in a local 
coordinate system (X 1 ) on M, we get 



Thus a critical curve satisfies 

d ( dL\ dL 



dt \dX l ) dX 



Fi , for i = 1 , . . . , dim M. 



We now assume that the manifold M is fibred over N, i.e. there is a bundle map n : M — >■ N. 
We denote the kernel of Tir by Vn C TM and is called the bundle of vertical tangent vectors. 
A cotangent vector a m G T^M can be restricted to V m ir, and such a restriction will be denoted 
by a v m : V m ^ — > K- Moreover, we denote the product bundle TN Xn M by TmN . Elements 
in this bundle arc denoted by (v n ,rn), where v n S T n N and m € M n = 7r _1 (n). The map 
(Ttt, t m ) ■ TM ->■ T M N; v m ^ (Tn{v m ),m) is denoted by p x . 

Definition 3 An intrinsically constrained Lagrangian system is a triple (jr : M — > N, L, F), with 
L a function on TmN and F the force term F : TM — ¥ T*M fibred over the identity on M . A 
curve m(t) G M is said to be a critical curve of the intrinsically constrained Lagrangian system 
if it is a critical curve of the associated Lagrangian system (M,p\L,F). 

It should be clear that the Lagrangian system [M,p\L,F) is singular since p\L is independent 
of the velocities in the fibre coordinates. We now focus on the Euler-Lagrange equations for 
the Lagrangian system (M,p\L, F) associated to the intrinsically constrained system (it : M — > 
N, L, F) (from now on, we denote L and p\L by the same symbol). For that purpose, let (x z , y a ) 
be a coordinate system on M that is adapted to the fibration, i.e. (x l ;i = 1, . . . , dim TV) are 
coordinates on N and (y a ; a = 1, . . . , dimM — dim N) are coordinates for the typical fibre of n. 
It is not hard to see that the Euler-Lagrange operator for an intrinsically constrained Lagrangian 
system, when restricted to vertical directions becomes 

8T 
£C{L)V ^ dy-^- 



A critical curve m(t) = (x 1 (i) , y a (t)) of the Lagrangian system on it : M —> N then satisfies the 
(possibly nonholonomic) condition 

dL 

-=-^(x,x,y) + F a (x,y,x,y) = 0, for a = 1, . . .,dimM - dimiV. (1) 

The appearance of this constraint on the set of critical curves justifies the previous definitions. 
We now specialize to three distinct cases for constrained Lagrangian systems. We focus on a local 
treatment. 



2.1 Gyroscopic forces 

We study here the specific case where the non-conservative force term F is of the form F(v m ) — 
—i Vm j3 m with j3 a closed 2-form on M. In this case we say that F is a gyroscopic force term 
associated to j3. As will become clear further on, Routh reduction of a Lagrangian system with 
non-abelian symmetry is in general a Lagrangian system subjected to gyroscopic forces. This 
type of forces therefore deserves special attention. Let v m — (x l ,y a , x 1 , y a ) in a coordinate system 
(x l ,y a ), then F assumes the following form 

F(v m ) = {faff + p ia y a ) dx l + (p ab y b ~ p. la x l ) dy a . 

We say that the gyroscopic force term is regular if /3 a (, is nondegenerate at all points in M, 
i.e. if is a nondegenerate 2-form when restricted to Vn. For gyroscopic forces, the intrinsic 
constraint §T§ becomes: 

dl 

^ + (/W-/3^)=0, 

and if /3 is regular, it should be clear that this constraint can be brought into the form 

r = - mv^rr" (^ - ^ ■ 

An interesting property of intrinsically constrained Lagrangian systems with a gyroscopic force 
term is that they admit a presymplectic formulation on the manifold TmN. To show this, we 
introduce the following maps: 



1. the Legendre map W X L : T M N -> T* 4 N, defined by 

L(v n + ew n ,m), 



{¥ 1 L(v n ,m),(w n ,m)) = — 

de 



e=0 



for (v n , m), (w n , 7Ti ) G TmN arbitrary; 

2. the energy El of the Lagrangian L is the function on TmN defined by 

E L = (¥iL(v n ,m),(v n ,m)) - L(v n ,m); 

3. the projections tti : TmN — > TN and TT2 '■ TmN — > M, and their analogues on the cotangent 
bundle: Wi : T M N -> T*N and 7f 2 : T^N -4 M. 

We denote the Liouville form on T*N by On and the associated symplectic 2-form lon = d0N- 
The following result is then easily proven in a local coordinate system. 



Proposition 1 Let m(t) denote a critical curve of the intrinsically constrained Lagrangian sys- 
tem (M, L, F) with gyroscopic force term associated to /3. Then j(t) = (h(t), m(t)) £ TmN, with 
n(t) = n(m(t)), satisfies the following presymplectic equation 

(i-V(t) ((7T1 o ¥xL)*uj n + ir* 2 ll) = -dE L ) | 7 . 

The converse also holds: any curve j(t) E TmN of the form ^(t) = (n(t),m(t)) with n{t) = 
7r(m(i)) and that solves the above presymplectic equation, projects to a critical curve m(t) = 
7r a (7(t)) of(M,L,F). 

If (i) F]Z : TmN — > T M N is a difFeomorphism and if (ii) /3 determines a regular gyroscopic force 
term, then the 2-form ((7F1 o¥iL)*ljn + k<20) is symplectic. In [12] the symplectic structure of 
a Routh reduced system is proven to be of this type and it is identified with the the Marsden- 
Weinstein reduced symplectic structure of a Lagrangian system with symmetry. Without the 
above regularity assumptions (i) and (ii) , the presymplectic system may only have solutions on a 
submanifold of TmN. A next step would be to find these submanifolds, however, this is not the 
scope of this paper and we refer the reader to the presymplectic constraint algorithm [6j [7] . 

We end this paragraph with a slightly more general formulation of the above proposition. In 
general, the intrinsically constrained Lagrangian systems obtained from Routh reduction will 
have a force term F = F\ + F2 that consists of two parts: a gyroscopic force term F\ associated 
to j3 and a force term F2 = T*ir o F o p 1; with F : TmN — > T*N. Again from local coordinate 
expressions, we conclude that critical curves of (n : M — ► N,L,F) solve the presymplectic 
equation (and vice versa) 



(i m ((¥1 o F!L)*u;jv + tt 2 */3) = -dE L + ttJ #) 



with rrlF : T M N -»■ T*{T M N) such that «f(»„,»n),y K , m )) = (F(v n ,m),TTn(V {Vn!m) )}, with 
(v n ,m) E T M N and V {Vn ^ m ) E T^ Vn>m )(T M N) arbitrary. 

2.2 Regular configuration constraints 

Throughout this paragraph we assume that the force term F is of the form F = T*ix o F o pi 
with F: T M N ->■ T*N as above, that is, F(v m ) = T; n ir(F(v n , m)), with v n = Tir(v m ). Then the 
constraint equation ([TJ reduces to 

dL/dy a (x,x,y) = 0. (2) 

Intrinsic constraints of this form are called configuration constraints and they are said to be 
regular if a map 7 : TN — > M, with n o 7 = r^v exists such that, in local coordinates, Eq. ([2} 
is equivalent to y a = 7 a (x,i). As we show here, in this case critical curves of the Lagrangian 
system on M are in a one to one correspondence to critical curves of a new Lagrangian system on 
N. The local existence of such a map 7 is guaranteed if the (vertical) Hessian of the Lagrangian 
is non-degenerate, i.e. if 

d 2 L 

dy a dy b 
Consider the Lagrangian system (N,L',F') on N, with 



det-— —r ^ 0. 



ll = L(v n ,j(v n )) and {F'{v n ),w n ) = {F(v n ^{v n )),w n ). 



Let m(t) = (x l (t) , y a (£)) be a curve in M. It is easily seen that the condition for m to be critical 
w.r.t (M — )■ N,L,F) is equivalent to n(t) = ir(m(t)) being a critical curve of (N,L',F'). More 
precisely, 



is equivalent to 



y a = l a (x,±). 



d (dV. .\ dL' .. _. . .. .. 



The correspondence is made explicit by the map 7: 7(n(t)) = m(t). In this sense the intrin- 
sically constrained Lagrangian system (M — > N,L,F) is 'equivalent' to the Lagrangian system 
(N,L',F')onN. 

The relation between Routh reduced systems and intrinsically constrained systems can now be 
made more explicit, and for this purpose we go ahead of some definitions and notations. We 
believe however that the reader familiar with geometric Routh reduction will benefit from the 
following. In general, a Routh reduced Lagrangian system will be an intrinsically constrained 
system on M — (Q/G^ x §) fibred over N = Q/G, where Q is a manifold on on which a Lie 
group G acts freely and properly. The bundle M is a product bundle, and the Routhian is a 
Lagrangian on TM independent of the velocities in the fibres of M. As will be made more precise 
in the following, this leads to constraints on the set of critical curves and due to the fact that M 
is a product bundle we will distinguish two constraints. The constraints arising from the bundle 
Q/Gfj, — > Q/G will always belong to the regular gyroscopic class. The constraints arising from the 
bundle q — > Q/G will be configuration constraints. If these configuration constraints are regular, 
then we retrieve the standard Routh reduction procedure. The irregular case will sometimes lead 
to a variational problem on N as is described in the next paragraph. In general, a Routh reduced 
system is associated to a presymplectic structure on TmN, and its solutions can be constructed 
by means of the presymplectic constraint algorithm. 

2.3 Linear constraints and Lagrangian multipliers 

Throughout this paragraph we describe a special case of configuration constraints for a Lagrangian 
system on M with the property that its critical curves correspond to the solutions of a variational 
problem on N. We assume that (i) w : M — > N is a linear fibration and (ii) the force term F is 
of the form F = T*ir o F o Ttt, with F :TN -)• T*N. 

We say that the intrinsically constrained system (M, L, F) is linear if L is of the form 

L(v n ,m) = L a (v n ) + (a(v n ),m). 

with a : TN —> M* , fibred over the identity and Lq a function on TN. The constraint equation ([1} 
for critical curves m(t) — (x l (t) , y a (£)) now becomes 

— (x,x,y) = a a (x,x) = 0. 

Assume that is a regular value of a, and consider the submanifold C of TN defined by v n s C if 
and only if a(v n ) = 0. The following proposition states that critical curves of (it : M — > N, L, F) 



coincide with critical curves of the variational problem on N with Lagrangian Lq constrained 
to the submanifold C in TN. Therefore, we can say that the initial Lagrangian system on M 
admits a variational interpretation on the base manifold N, independent of the total space M. 
When working in a coordinate system, one might interpret the fibre coordinates y a as Lagrangian 
multipliers for the constrained variational problem on the base manifold. 

Definition 4 A critical curve p(t) of a Lagrangian system (P,L',F') on a manifold P con- 
strained to a submanifold C of TP is a curve satisfying 

8 J L'(p(t))dt + J(F'(p(t)),Sp(t))dt = 0, 

for any variation 8p with fixed endpoints originating from a deformation p e (t) of p(t) which is 
entirely contained in C , i.e. p e (t) £ C . 

Proposition 2 Given a linear intrinsically constrained Lagrangian system, (n : M — > N, L = 
Lq + (a, m), F) on M . Then the critical curves of (tt : M — > N, L, F) on M project under ir onto 
critical curves of the Lagrangian system (N,Lq,F) on N constrained to C . 

Proof. The proposition is easily proven by noting that critical curves of (ir : M — > N, L, F) 
satisfy 



S / L(m{t))dt+ / (F(h{t)),5n(t))dt = 

for arbitrary variations with vanishing endpoints in M and project onto curves in C. Among 
all these variations, we can now consider variations originating from deformations m e (t) of m(t) 
that project onto deformations n e (t) of n(t) that are contained in C. In that case the first term 
in the latter equation is precisely 



d_ 
~dc 



L(rhJt))dt = — 
de 



L Q (h e (t))dt, 
and we thus find that n(t) is a critical curve of (N, Lq, F) constrained to C. 

2.4 Connections 

We now return to the general setting of intrinsically constrained Lagrangian systems and conclude 
this section with a proposition which describes how to split up the Euler-Lagrange equations 
making use of a chosen connection on the fibration 7r : M — > TV. A connection on M is a direct sum 
decomposition of TM, i.e. a distribution H on M such that TM = H © Vw, and the associated 
map h : TmN — > H is called the horizontal lift. Using this connection we are now able to say that 
a variation 8m(t) of m(t) is horizontal if 5m(t) £ H, or if equivalently 5m(t) = h(m(t) , 5n(t)) , 
where Sn(t) is the corresponding variation of n(t) = 7r(m(i)). We denote the horizontal and 
vertical components of an arbitrary tangent vector w by w h and w v respectively. The restriction 
of a cotangent vector a to the horizontal distribution is denoted by a h . 

Proposition 3 The Euler-Lagrange operator of an arbitrary Lagrangian system (M, L, F) on M 
is completely determined if one only considers horizontal and vertical variations, i.e. if m(t) is 
a critical curve of (M, L, F) then it satisfies: 

6 / L{m(t))dt + / (F(m(t)), 5m(t))dt = with 8m v (t) = 0, 
S ( L(rh(t))dt + f(F(rh(t)), 8m(t))dt = with Sm h (t) = 0. 



The first and second equations are equivalent to, respectively, 

EC{m{t)) h + F{m{t)) h = 0e T*N = H m , 
£C(m(t)) v + F(m(t)) v = G V£>r. 

Assume an adapted coordinate chart is chosen and that the horizontal distribution is spanned by 
the following set of vector fields 

dx> LAX ' V) dy°' 

with r° the connection coefficients, then the horizontal and vertical parts of the Euler-Lagrange 
equations take the following form 

w ~ di {d¥j +Fi ~ [dr + Fa ) rf = °' 

dy a 

In fact, by substituting the second in the first equation, we may conclude that these equations 
are independent of the connection coefficients. 



3 Lagrangian systems with symmetry 

We consider a Lagrangian system (Q,L,F) on a manifold Q and we assume that a Lie group 
G freely and properly acts on Q. This ensures that Q/G has a manifold structure and that 
7r : Q — > Q/G is a principal fibre bundle. The (right) action is denoted by \I/ : G x Q — > Q : 
(.9i <z) ^ ^(li 9) = 19- The associated maps a q : G — ?> Q and R g : Q — ?> Q are defined by 

qg = o- q {g) = R g {q). 

For the sake of simplicity, the infinitesimal action induced by a will be also denoted by a: 
a q(0 = Ta q (t;), with £ e q a Lie-algebra element of G. The vector field £q on Q is defined by 
£q(?) = a q(0- We can consider the lifted action of G on TQ, and using the previous definition 
on TQ, we shall write £tq to denote the symmetry vector field associated with the Lie algebra 
element £ G g. Throughout this section we shall extensively rely on the choice of a principal 
connection on Q. This is a g- valued one-form u on Q satisfying w(£q) = £ and which is equivariant 
in the sense that, for arbitrary g G G, 

R*u> = Ad g -i ■ u, 

where Ad denotes the adjoint action on g. We assume that the notions of horizontal lift, associ- 
ated bundles and covariant derivatives, etc are well-known (see [ID]). One typically defines the 
momentum map associated with the Lagrangian L in the following way. 

Definition 5 The momentum map Jl '■ TQ — y g* is the fibre derivative of L restricted to vertical 
directions: 

L(v q + eZ Q (q)). 

£=o 



Jl(vM) = j e 
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It is well-known that if the Lagrangian is invariant under the (prolonged) action of G on TQ, 
then (i) Jl is an equivariant function, i.e. R*Jl = Ad* • Jl and (ii) Jl is a first integral of the 
critical curves of the conservative Lagrangian system (Q,L,F = 0). 

Definition 6 A force term F : TQ — > T*Q is said to be invariant under the action of G (or 
simply G-invariant) if the following two conditions are satisfied: 

L {F{v q g),w q g) = (F(v q ),w q ) 

for v q ,w q € T q Q and £ £ g arbitrary. 

A Lagrangian system (Q,L,F) is G-invariant if the Lagrangian L is an invariant function under 
the action of G and if the force term F is G-invariant. 

A straightforward consequence of the definition of a G-invariant force term F is that F is reducible 
to a map / between the quotient spaces TQ/G and T*(Q/G). Moreover, the second condition 
in definition [5] guarantees that the momentum map Jl is conserved along critical curves of a 
G-invariant Lagrangian system (Q,L,F) (it is sufficient, but not necessary). In the following 
proposition we study the necessary conditions for a Lagrangian system (Q,L,F) such that Jl 
is conserved along its critical curves. This generality will be of importance further on when 
considering the Routhian and its momentum map, both associated with an invariant Lagrangian. 



Proposition 4 The momentum map for a (not necessarily G-invariant) Lagrangian system 
(Q,L,F) is conserved if and only if the force term is such that (cIL,£tq) = — (F,£q) fo r a ^ 
£ G g along the critical curves. If this holds we say that the Lagrangian system (Q,L,F) has 
conserved momentum Jl w.r.t. the action of G. 

Proof. Assume that q(t) is a critical curve of the Lagrangian system (L,F). We now compute 
the derivative of the momentum map along the trajectory and proof that it is equal to zero. Let 
(q l ) be a local coordinate system on Q and let £ G g be chosen arbitrarily, then 



d_ 
dt 



(Mmm) 



d_ 
dt 



d_ 
de 



-={) 



L(q(t) + et Q (q(t))) = - 



dL\ 



dt \dq l )^ Q 



Mi 

dq 1 ^ 



FiC Q 
$TQ(L)(q(t)) 4 



dLd% 

— r — -Q 
dq l dq> 

, OL dfy 

I — —q 

dq % dq J 

(F(q(t))^ Q (q(t)))=0. 



dLi 

dq^ Q 



Note that, if L is invariant then the momentum map is conserved along critical curves if F 
satisfies condition (2) from definition [6] 

The first step in our approach to Routh reduction of a G-invariant Lagrangian system (Q, L, F) 
consists of defining an equivalent Lagrangian system on Q whose Lagrangian is only invariant 
under the the action of a subgroup of G although it has conserved momentum w.r.t. the action 
of G. This observation should justify the above definitions. 
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3.1 General constructions on the quotient spaces 

Assume throughout the remaining of this paper that a regular value fi 6 g* for the momentum 
map of a G-invariant Lagrangian system (Q, L, F) is chosen. We denote the group G M as the 
isotropy subgroup of G of ^ under the coadjoint action of Gong*. 

This section describes the structure of the manifold TQ/G ^ on which the Routhian lives and 
how variations of a curve in Q behave when projected onto the quotient TQ/G^. Many of the 
properties explained below can be found in more detail in [2]. We start with fixing notations 
and defining several projections and bundles appearing in the chosen framework. Recall that the 
bundle Q/G^ — > Q/G is an associated bundle of the principal fibre bundle Q — ¥ Q/G, whose 
typical fibre is the orbit O^ of fj, under the adjoint action on g*. 

Definition 7 The projections onto the different quotient spaces are denoted by 

n:Q^ Q/G, n^-.Q^ Q/G,,, 
Tfp : Q/G^ -»• Q/G, n:~ 9 :=(Qx g)/G -)■ Q/G. 

A point in Q, Q/G '„ and Q/G is typically denoted by q, y and x respectively. Similarly, tangent 
vectors in T q Q, T y (Q/G^) and T X (Q/G) are denoted by v q , v y and v x respectively. Further, 
elements in the bundle g — > Q/G associated to the Lie-algebra g are denoted by £ = [q, £]q, with 
£ e g, q e Q. 

If a principal connection u) on 7r : Q — > Q/G is chosen, then the horizontal lift of a tangent vector 
in T{Q/G) to the point q 6 tt^ 1 (x) is denoted by (v x ) q . The horizontal lift of a curve x(t) from 
Q/G to Q is denoted by qh(t) 

According to our previous notation conventions, an element in the product bundle T(Q/G) x 
Q/Gy, xg^ Q/G is denoted by (v x ,y,£), where v x € T X (Q/G),y e Q/G^ and | <E g satisfy 

*»(v) = m) = xeQ/G. 

Proposition 5 The manifold TQ/G^ is isomorphic to T(Q/G) X Q/G^ x g. The isomorphism 
can be made explicit with the choice of a principal connection ui on Q. 

Proof. Assume that a connection u) on Q is chosen. Let [i\j]g, £ TQ/G^ be the orbit of an 
arbitrary tangent vector v q in TQ. We now define a map 

^ ■ TQ/G^ -*■ T(Q/G) x Q/G^ x g, (3) 

together with its inverse tp u as 

i(K]cJ = {T-K(v q ),7r^(q),[q,u}(q)(v q )] G ), 

ipcu(v x ,y,0 = [(v x )g +(Tq(Q]Gn, withy= [q] Gfl and f = [q,£] G . 

It is standard to show that these maps are smooth and do not depend on the choice of an element 
in the orbits under consideration. 

Using similar arguments as before, one can show that TQ/G = T(Q/G) x g. Taking into account 
these isomorphisms, the natural projection A M : TQ/G^ — > TQ/G simply reads 

*»(vx,v,£) = (y x ,£). 
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Definition 8 An invariant Lagrangian function L on TQ reduces to a function I on TQ/G, 
defined by 

l(v x ,i) = L((v x f q + a q (0), with q e tt-^x), [q,£] G = f. 

The function X*l, defined on TQjG ^ is denoted by the same symbol. 

The associated momentum map Jl '■ TQ -4 g* , reduces to a bundle map ji : TQ/G — > g*, fibred 
over the identity on Q/G and defined by 

{jl( v x,0,fj) = {Jl{{v x )\ + o- q (£,)), n), withq € ^(x), [q, £] G = £,[q,rj\ G = fj. 

The reduced map ji can alternatively be defined from I directly. Let FA : T(Q/G)xg — >• T(Q/G) x 
g* be the fiber derivative of I with respect to the (linear) factor g. 

Proposition 6 The maps ji and I are related by: 

Fil(v x ,i),(v x ,fj)) = {ji(v x ,t),ij). 

Proof. Let (v x ,£), ( v x,v) £ T(Q/G) x g be arbitrary. Fix a point q 6 Q projecting to x such 
that [q,£] G = I and [q,rf\ G = fj; then 



FiHv x ,i),{v x ,ti)) := ^ 



l(vx,£+£fj) = -f 

£=0 ae 



L((v x ) h q +<j q (Z + en)) 
i 

The element (i £ g* defines a map ft : Q/G^ — > g* such that 

</*(»), I) = (A*, 0. with [g]c„ = y and [«, £] G = £. (4) 

Due to the equivariance of Jl, the level set J£~ (/x) is G^ invariant and reduces to a subset of 
TQ/G^. 

Proposition 7 The quotient manifold J^ (^)/G /J is identified with the subset of points (v x ,y,£) 
in TQ/G,, satisfying ji(v x ,£) = ft(y). 

In classical Routhian reduction, an essential condition for reducibility of an invariant Lagrangian 
system is that this submanifold j£" (n)/G^ is precisely T(Q/G) x Q/G^. This situation is the 
subject of the next session. 

3.2 The regular case 

Definition 9 An invariant Lagrangian L is said to be G-regular if the map Wd is a diffeomor- 
phism between T(Q/G) x g and T(Q/G) x g* . 

We shall write the inverse of this map as ki : T(Q/G) x g* -4 T(Q/G) x g. 

Proposition 8 If L is G-invariant and G-regular then the mapping 

JZ\lS)/G» -4 T(Q/G)xQ/G^ 
K]g„ >-► (TTT(v q ),[q] G ^ 

is a diffeomorphism. 
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Proof. We provide the inverse of the mapping of the statement. Given a point (v x ,y) £ 
T{Q/G) x Q/Gp, and let £ £ Q be such that 

(«*>£) = ni(v x ,fi(y)). 

Fix a representant g for y, i.e. y = [?]g„, and let £ be such that £ = [<7, £]g- Then the point 
w g = (( w a;)g +£q(<z)) G ^Q belongs to the set J z ^ 1 (/i). Indeed, with ry = [q,??]g arbitrary, 

Ji((vx), +v q {Q){v) = (ji(vx,£),fj) = (fj.(y),v) = (ji,n). 

It only remains to show that this definition is independent of the chosen representative q at y. 
Due to the equivariance of Jl this is standard. Finally, with (v x ,y) we can now define an element 
in J^ 1 (fi) J G fj, determined as the orbit of v q . The mapping from T(Q/G) x Q/G^ to J£" 1 (/i) 
obtained in this way is the inverse of the map defined in the statement of the proposition. 

Remark 1 Note that the mapping given in Proposition [S| is precisely the first two components 
of the diffeomorphism {3J) in Proposition [5l In other words, the G-regularity shows that the 
submanifold JJ) (/k)/G m in T(Q/G) x Q/G^ x g can be given as the set of points {v x ,y,£) with 
{vxilj) £ T{Q/G) x Q/G^ arbitrary and with £ such that (v x ,£) = ki(v x , A(y)). 

With this in mind, it follows that the converse of Proposition^ is also valid: if for any p £ g* 
the map JJ) (^t)/G A , —¥ T(Q/G) x Q/G^ is a diffeomorphism, then L is G-regular. This shows 
that our definition of G-regular Lagrangians is independent of the connection u). 

Remark 2 Let p be a G -invariant Riemannian metric on Q and V £ C°°(Q) a G -invariant 
potential energy. The Lagrangian L{v q ) — ^p(v q ,v q ) — V(q) (kinetic minus potential) is G-regular 
(see, for example, U5f). 

Note however that a Lagrangian can be both hyperregular (i.e., its Legendre transformation is a 
diffeomorphism) and G-invariant but not G-regular. Recall the example from the introduction 
L = (q 1 ) 2 + q l q 2 - V{q l ) with V £ G°°(R), Q = M 2 and G = R acting by translations in 
the q 2 coordinate. Fix a connection, say u> = dq 2 , and let £ = {q 1 ,£,),fi — {q l ,ii) be arbitrary in 
= Q/GxR. Then the reduced Lagrangian I : T(Q/G) xq ->• R is l{q l , q l ,£) = (q^f+q^-Vlq 1 ) 
and (F/;(g 1 , q , £), (q ,q ,ff)) = q 1 ^. Thus ¥1? is not a diffeomorphism. 

It is standard in Routh reduction to introduce the following definition, although due to the fact 
that L may not be of type 'kinetic minus potential' the definition is slightly more general than 
the usual definition of the locked inertia tensor. 

Definition 10 Given a Lagrangian L : TQ — > R and a point v q £ TQ, the mapping T v : g — >• q* 
defined by 



z»,(0fa) 



d 

Te 
d_ 



d_ 



L(v q + e£ Q (q) + Tr, Q (q)) 

T = 



e=0 

J L (v q + r VQ (q))(0 



for £,?7 £ q arbitrary, is called the locked intertia tensor at v q . 

We say that L is locally G-regular if the locked- inertial tensor I v is invertible for any v q £ TQ. 
Local G-regular Lagrangians only guarantee that J£ 1 (fj.) is locally diffeomorphic to T(Q/G) x 
Q/G^. In any case, local G-regularity presents the advantage of being a property easy to check 
(besides the applications of the locked inertia tensor in other situations as stability, cf. [S]). We 
refer to [3] for a description of Routh reduction for locally G-regular Lagrangians. 
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3.3 The reduced form induced by the connection 

We now return to the general case and continue denning additional objects living on the quotient 
spaces under consideration. These objects will appear as force terms when studying Routhian 
reduced systems. Recall that we fixed a principal connection uj on Q. The one-form on Q obtained 
by pairing /i with the g- valued one- form w is denoted by w M — (/-t, w). The Lie-algebra g M of the 
isotropy subgroup G M is a subalgebra of g which can alternatively be determined as the set of 
elements (eg satisfying adj/i = 0. 

Proposition 9 The two-form dio^ — (fi, dui) on Q is protectable to a two-form on Q/G^, denoted 
by^. 

The invariant force form F : TQ — > T*Q reduces to a map: f : TQ/G — > T*{Q/G). 

Proof. It is easily seen that dcu^ is invariant under G^ and annihilates vertical tangent vectors 
in Q — > Q/G^. The latter follows from the equivariance of uj: let £ be an arbitrary element in 
Qn, then 

= (/j,, — adjw) = — (adt/x, w) =0 . 

The second part of the proposition is a straightforward consequence of Definition [5] of an invariant 
force term. Next we study the bundle W M : Q/G^ — > Q/G and in particular the set of tangent 
vectors vertical to 7f^. 

First note that both it : Q — >• Q/G and n^ : Q — > Q/G^ are principal fibres bundles (from 
now on PFB), with structure group G and G^ respectively. The adjoint bundle g is the bundle 
associated to g w.r.t the PFB x. Similarly we can consider the adjoint bundle g^ associated to g^ 
w.r.t. the PFB n^. It is not hard to see that g M is a subbundle of the fibred product Q/G^ x g: 
let [<?,£] g G 0m be arbitrary, then we define an element in Q/G^ x g by fixing a representative 
(q,0 S Q x g^ and we let (y = x^{q) — y,£ = [<?, £]g) determine the corresponding element in 
Q/Gfj, x g. The Lie algebra bracket [■, ■] on g carries over to a Lie algebra structure on the fibres 
of g. The bracket is denoted with the same symbol: [£,77] = [q, [£,t?]]g, with [<?,£]g = £ and 
[q, t]]g = fj arbitrary. For notational convenience we introduce the following shorthand notations: 
Tf^fl = Q/Gfj, x g, and the dual bundle by 7r*g* = Q/G^ x g*. Since g M is a subbundle of 7r*g 
we may consider the quotient bundle: x*q/q^. Finally, note that ft : Q/G^ — > g* can be seen as 
section of Q/G^ xg*-> Q/G p. 

Definition 11 Associated with the map fi : Q/G^ — > g* we define a section ad*/i of /\ tt*q* — > 
Q/G^ in the following way 

&d*fi(y)((y, 0, (y, fj)) = (fi(y), [£, fj}), 
where £,?y £ Qx,^p,{y) = X arbitrary (recall that 7f*g = Q/G^ x g). 

We will use the following shorthand notations: ad*/2(y)(£, fj) or ad|/i(y) = i?(ad*/2(y)). The two- 
form ad*/2 is clearly antisymmetric and its kernel consists precisely of elements in g^. Therefore, 
without introducing a new symbol, we can regard ad*/2 as a section of /\ (7f^g/g M )* — > Q/G^ 
which is nondegenerate as a map 7f*g/g p — > (7f*g/g M )*. Recall that by choosing a point q e 
7r u 1 (y)i * ne fibre of the subbundle g M at the point y is isomorphic to g M , the Lie algebra of the 
isotropy subgroup G^. This observation brings us to the next proposition. 
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Proposition 10 The bundle Vir^ — > Q/G ^ of vertical tangent vectors to the fibration tt^ : 
Q/Gfj, — > Q/G is isomorphic to TT*^g/g^ — > Q/G M . 

Proof. Wc define a surjective linear map from ~n*g to Vtt^, in the following way. Fix a point 
(y,£) £ 7f*g such that y — [q\c tl and ^(j/) = x - Let £ be a representant of £ at the point q, 
i.e. £ = [g, £]g, and consider the vertical tangent vector in VyW^ obtained by T7r M (<7q(£)). It is 
clear that this map is well defined (i.e. independent of the choice of the point q in [q]c and 
that it is onto. Furthermore, the kernel of this map consists of the elements in g^ at the point 
y. For later purposes it is important to see that the tangent vectors to the vertical curves 
[95(*)]g m e 7f~ 1 (a;) through y with g(0) = e and £ = g(0) £ g are mapped onto (y, [q, £] G ) + M at 
the point y. We conclude with recalling the definition of a connection on an arbitrary associated 
bundle corresponding to the fixed principal connection to on Q. Let E denote a manifold on 
which G acts from the left. 

Definition 12 The connection or horizontal distribution on the associated bundle Q Xg E cor- 
responding with the principal connection to is defined as the set of tangent vectors to curves of 
the form [q(t),f]c with q(t) £ Q horizontal and f £ E arbitrary. In the case that F is a linear 
space on which G acts linearly, we define the covariant derivative D as 

— [ q (t),e(t)} G = [q(t), (u} g(t) (q(t))) ■ e(t)] G + [q(t), e(t)] G , 

It is understood that in the above definition the action of a Lie-algebra element to q (q) on the 
element e £ E is denote by •. 

We will use the connection in associated bundles in the following two cases: 

Firstly, the principal connection to determines a connection in Q/G^ — » Q/G, the bundle asso- 
ciated to the orbit space C M . The tangent space to Q/G^ can be written as a direct sum of the 
horizontal and vertical distribution: 

T(Q/G„) S 7f;T(Q/G) ® Q/G „ %m»- 

This in turn implies that the two-form f3^ (cf . Proposition [9] above) can be decomposed into a 
horizontal-horizontal and vertical- vertical part (the horizontal-vertical part vanishes) 

/^ = (fF,-ad*/i), (5) 

where 

0" : Q/G^ A 2 T*(Q/G);y^ {£&),«>, 

fi being the g-valued curvature 2-form on Q/G. Such a decomposition is a straightforward 
consequence of the structure equation fi = duo + [lo, lo] and the fact that in this specific case the 
horizontal distribution on Q/G^ is the projection of the horizontal distribution on Q. 

Secondly, if E = g the covariant derivative on the adjoint bundle g equals, for any q(t) £ Q and 
£(£) £ g (with [•, •] the Lie bracket on g): 

2.[q(t),at)h = [<?(*), hKt)(?(*)U(t)]] G + WU(t)]G. 
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3.4 Variations 

The notions introduced here are described in more detail in [2J. For the sake of completeness 
and because we work in a slightly different situation we nevertheless briefly introduce the notions 
of vertical and horizontal variations in the context of Lagrangian systems with symmetry. The 
main point is to consider a variation of a curve q(t) in the total space Q, and subsequently study 
the projection of this variation to the quotient spaces Q/G^ and Q/G. 

For that purpose, we fix a curve in q(t) in Q and its lift q(t) in TQ. Following the notations from 
the previous section, we shall write the projection of q(t) to TQ/G^ as {x(t),y(t),£(t}), where 
x{t) = T q(t) Tr(q(t)), y(t) = 7r p («(t)) and £(t) = [q(t),u q(t) (q(t))} G . 

Definition 13 A deformation of the curve q(t) is smooth function q(t,e) — q e (t) such that 
q(t,0) — q(t). The corresponding variation is defined by 



Sq(t)= dq{t ' e) 



ck 



£ T q (t)Q 



— 



1. A vertical deformation is a deformation that can be written as q(t,e) — q(t)g e (t), where the 
family of curves g e (t) £ G satisfies go(t) = e, with e the unity of G. 

2. A variation dq(t) such that Tir q ( t )(5q(t)) — for all t is said to be a vertical variation. 

3. A variation Sq(t) such that u)(q(t))(8q(t)) — for all t is said to be a horizontal variation. 

Obviously, the variation induced by a vertical deformation is vertical and it is not hard to see 
that any vertical variation can be given in that way. 

Given an arbitrary vertical deformation q(t, e) = q{t)g f {t), then the projection of (q(t, e), q(t, e)) 
to the quotient space TQ/G^ is written as {x{t), y(t, e),£(£, e)), where according to the previous 
definitions: 

y(t,e) = [q(t,e)] Gll 

|(t,c) = [q(t)g e (t),u(q(t,e))(q(t,e))] G 

= [q(t)g e (t),u(q(t)g e (t)) (TR gAt) (q(t) + <J q {g e g7 1 )))] G 

= [q(t),u(q(t))(q(t)) + 9 e(t)g- 1 (t)] G - |(t) + [q(t), g e (t)g-\t)} G 

The typical structure equations of reduced vertical variations is obtained in the following 
way. Let d/de e=0 5e(0 = 8g(t) s g and fj(t) = [q(t),Sg(t)]c S fl- Using the definition of the 
covariant derivative on the associated bundle § — >• Q/G, we may write 

^fj(t) = [«(*), [u J ( q (t))(q(t)),Sg(t)]] G + [q(t) : 6g(t)} G . 

Finally, by computing the derivative to the curve e — s- £(£, e) we find the structure equation [2]: 



9x(t) 3 6£(t) := — 



£(t,e) = [q(t),5g(t)} G = -^ry(i) - [£(t),ij(t)]. 



It now remains to check the variation in the variable y. The curve e i-> y(t,e) = [q(t,e)] G is 
contained in a fibre of the bundle Tf^ixit)) and consequently, the tangent vector to this curve is 
in VW^. In the language of Proposition [TU1 the variation 8y in y then corresponds to (y, ff) + q^. 
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Next, we describe the structure of reduced horizontal variations. Let q(t, e) denote a defor- 
mation, with corresponding horizontal variation 8q(t). Again, the projection of the deforma- 
tion q(t, e) determines a deformation (x(t, e),y(t, e),£(£, e)) of (x(t),y(t),t;(t)) in TQ/G^ (with 
x(i, e) = w(q(t, e))). By definition we have that £(i, e) = [g(i, e),w(q(t, e))(q(t, e))]c- Using the 
fact that 5q is horizontal, i.e. Sq = 5x^, the equality 

-£ «(«(t, C ))(g(t,c))=n(«(t))(d«(t),g(t)) 

ae 

holds and we are able to write the horizontal and vertical parts of the tangent vector 5£(t) £ ^e(t)9 
as follows: 



*»■''■ = Is 



«MI = («0,[9((),S)(«(())(«(),9(())]o 

e=0 



= [l;(t),n(x(t))(5x(t),x(t)) 
Z(t,e)\ = {i,6x(t)), 



d 



5 ^ k = • rf e 

where we used the decomposition of TVmfl — tt*T(Q/G) © Vn into its horizontal and vertical 

subspace (due to the linear structure, Vtt is identical to g x 5). It is important to see that <5£ is 
not necessarily horizontal w.r.t the connection on the bundle g — > Q/G, despite the fact that it 
originates from a horizontal variation of q(t). 

The deformation y(t, e) has a corresponding variation Sy(t) which is a tangent vector along y(t). 
It is horizontal in terms of the connection associated to w on the associated bundle Q/G^ — > Q/G. 

3.5 Lagrange-Poincare reduction 

From the above structural equations for projected vertical and horizontal variations of a curve 
q(t), we are now able to deduce the Lagrange-Poincare equations on TQ/G = T(Q/G) x jj for 
an invariant Lagrangian system (Q,L,F) on Q. This follows easily from the following identity: 
given an arbitrary deformation of q(t), then 



L(q{t,e))dt= / l(x(t,e),£(t,e))dt. 
1 Ji 

Now, if q(t) is critical for L, then for any vertical variation with vanishing endpoints we write, 
with a slight abuse of notations 

= jjmt,e)) = jjmtu{t^))dt 

\-§- t ¥ i l - &d l ¥ i l 'V( t ))dt, 

for all fj(t) G g with vanishing endpoints. On the other hand if we consider horizontal variations, 
then 

= ^-[l(x(t,e),£(t,e))dt+[(f(x(t),$(t)),6x(t))dt 

((£C{l) h + f,6x) + (¥ i l(x,0,Cl(x)(8x,x)))) dt. (6) 
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Recall that I can be regarded as an intrinsically constrained Lagrangian on M = § (and N — 
Q/G), and hence ££{l) in the previous equation is the Euler-Lagrange operator for I regarded as 
a function on Tg. The reduced equations now read 

— ¥A = -ad*~F,~Z, and 
Dt ? « ? 

££{l) h = (F|UiO)-/. 

It should be clear that these equations are not the Euler-Lagrange equations for the intrinsically 
constrained Lagrangian system (jj — > Q/G, I, f) on the fibred manifold jjj — > Q/G. The introduc- 
tion of the Routhian will avoid this obstacle, paying the price that additional non-conservative 
forces have to be taken into account. 



4 The Routhian reduction scheme 

Consider a critical curve q(t),t £ [a, b], of the invariant Lagrangian system (Q,L,F), that is 

= 6 J L(q(t))dt = - J(F(q(t)),5q(t))dt, 

with Sq(t) an arbitrary variation of q(t). By considering variations with fixed endpoints, we 
obtain the Euler-Lagrange equations for critical curves 

£C(L)(q(t)) = -F(q(t)), Vt e I. 

We assume throughout this section that we fixed a principal connection uj and a regular momen- 
tum value fi G 0* of Jl- 

Definition 14 The Routhian is the function on TQ defined by 

R»(v q )=L(v q )-(n,uj(q)(v q )). 

The Routhian depends on the choice of the connection and the momentum fi. Although the 
Routhian is not in general invariant under the action of G, we do have the property that the 
momentum map of the Routhian w.r.t the action of G will be conserved (cf. propositions]) . We 
thus study the behavior of the Routhian under the action of G. 

Proposition 11 The Routhian is invariant under the action of G M . It transforms under the 
infinitesimal action of G as: 

Ztq (i? M ) (v q ) - (ad^, «(«)(«,)), V£ e g. 

Proof. This is straightforward from the definition of i? M and the equivariance of u>. 

Proposition 12 

1. A critical curve of the invariant Lagrangian system (Q,L,F) is a critical curve of the 
Lagrangian system {Q^R^^F + G M ) and vice versa, where G M is the gyroscopic force term 
(see '<\2.1\ above) associated to diu^ , that is 

G^{v q ) = -i Vq d^ 1 = -(n,i Vq du>), v q <E T q Q. 
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2. The Lagrangian system (Q, R 11 , F + G^) has conserved momentum Jrh w.r.t. the action of 
the entire symmetry group G. 

3. The momentum maps of L and R^ , both considered with respect to the full action of G, are 
related by Jrm = Jl — £i- Hence q(t) is critical curve of the invariant Lagrangian system 
(Q,L,F) with momentum Jl — fJ- if and only q(t) is a critical curve of Lagrangian system 
(Q, i? M , F + G M ) with momentum J^ = 0. 

Proof. The first statement is proven by considering the variation of i? M along a critical curve 
q(t) of (Q, L, F): given an arbitrary variation 5q(t) of q(t) in Q then 



S / R"(q(t))dt = S / L(q(t))dt-S / (w"(ff(t)),?(t))<ft 

= 5 J L(q(t))dt + J \i m d^(q(t)),5q(t))dt - (u»(q(t)), 6q(t)) b a 

= - J(F(q(tj), Sq(t))dt - J(GP(q(t))> Sq(t))dt 

-(u»(q(t)),Sq) b a + (¥L(q(t)),Sq(t))\ b a ■ (7) 

This easily proves the first statement in the proposition (the other direction follows by reversing 
the arguments). 

For the second statement, according to Proposition 21 it is sufficient to show that £,TQ(R >i )(v q ) = 
— (F + G M , £q )(q). The contraction of the force term F + G p with vertical directions is precisely 

= (m> [<*>(?) (« g ),£]) 
= -(ad£/i,w(g)(w g )). 

Together with the previous proposition this shows the second statement. 

Finally, the third statement is proven by computing the momentum map w.r.t action of G of the 
Lagrangian system (Q, R^, F + G M ): 



(jR»(Vq),0 = -r 



e=0 



L(v q + ea q (0) - (n,oj{q)(v q + ea q (£))) 



= (JL(Vq)-H,0- 

We are now ready to state the first theorem of this paper. First we need to fix additional 
notations. The Routhian is reducible to a function 

Vf : TQ/G^ 2 T{Q/G) x Q/G^ x § ->■ K. 

Recall that the invariant force term F reduces to a map 

f:TQ/G^T*{Q/G) 

and that the gyroscopic force term G^ reduces to a gyroscopic force term 

C:T(Q/G^)^T*(Q/G^) 
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on Q/G M associated to /3 M (see Proposition ^. With these notations and those in Sj5]m mind, we 
now consider the bundle M — > Q/G, where 

M := Q/Gn x § -> Q/G. 

The triple (M — » Q/G, 7^/*, / + C) w iU determine an intrinsically constrained Lagrangian system 
with 

T M N = T(Q/G) x Q/G^ x . 

To reduce the notational complexity, it is understood that / and £ M are pull-backed to the 
appropriate bundles in order to fit the definition of an intrinsically constrained Lagrangian system. 

Due to the product structure of the bundle M — Q/G^ x §, we can restrict the Euler-Lagrange 
operator ££(7\L M ) to tangent vectors vertical toM-ig or to M — > Q/G^. The latter is denoted 
by dylZ^ and the first is precisely F;7?. M . The two bundles Q/G^ and g are equipped with a 
connection determined by u>, and it is standard to show that these determine a unique connection 
on M for which the horizontal distribution projects onto the horizontal distribution of Q/G^ and 
0. It should be clear what is meant when we write £C(Jl^) h and £C(1l' M ) v = (d y 1l^ ,V /R,^) . 

Theorem 1 (i) A critical curve of the invariant Lagrangian system (Q, L, F) with momen- 
tum n projects onto a critical curve of the intrinsically constrained Lagrangian system (M — > 
Q/G, TV 1 , f + C M )- (ii) Conversely, if a critical curve of the intrinsically constrained Lagrangian 
system (M — > Q/G,1Z^, f + £ M ) is the projection of a lifted curve q(t) in TQ, then q(t) is a 
critical curve of (Q, L, F) with momentum \x. 

Proof. The proof of this theorem is obtained by first considering the Lagrange-Poincare re- 
duction of the Lagrangian system {Q, R^ , F + G M ) w.r.t the action of G M , and consequently by 
showing that the reduced equations for critical curves with zero momentum coincide with the 
critical curves for the intrinsically constrained Lagrangian system (At — > Q/G, W 1 , f + £ M ). 

Let q(t) : L — > Q denote a critical curve of (Q, L, F) with momentum \x. From proposition [T2l we 
know that q(t) is a critical curve of the Lagrangian system (Q, R^, F + G^) with momentum 0. 

The projection of q(t) onto M is denoted by (y(t), £(£)). Recall that the vertical part of y(t) equals 
(y W; £(*)) + 0m (^ ms i s the reduced version of the fact that the curve in TQ is the complete lift 
of a curve in Q, see infra). By definition of TU 1 , we may write i?^(g(i)) = lZ^(x(t), j/(£), £(£)), 
with x(t) = n(q(t)). We now show (i) by studying the reduction of the variational equation for 
critical curves of (Q, R^, F + G^) with zero momentum (similar to section ^3.5l for the Lagrangian 
system (Q,L,F)). Consider an arbitrary variation 6q(t) of q(t) and the corresponding reduced 
variation (5y(t),S^(t)) of (y(£), £(£)), then (for notational convenience we omit the explicit time 
dependence) 



5 / Rf(q)dt + ((F + G»)(q),5q)dt 

= 8 j 'n»(x,y,£)dt + J \(f(x,i),5x) + (e(y),Sy))dt. (8) 

Vertical Variations: We apply the Lagrange-Poincare reduction scheme and we assume Sq(t) to 
be a vertical variation. From the previous section, we know that the reduced variations Sy(t) of 
y(t) and <$£(£) of £(£) both satisfy the following equalities: 

s£(t) = §- t m-m,m} 

6y(t) = (v,fj(t))+Bn 
Sx(t) = 
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where fj(t) can be chosen arbitrarily (with or without vanishing endpoints). From Eq. ([5)1. we 
conclude that 

for reduced vertical variations of (y(t),^(t)). Moreover, since the variations are vertical w.r.t the 
fibration M — > Q/G, the first term is 



S / Wfay&dt = / (d v 1Z",5y) + (F^ l ,50dt. 

From the fact that F^T?.^ = ¥A — ji = 0, the variational equation holds for arbitrary vertical 

variations (5y, <5£) of (y, £), and not only those satisfying 5£ — j^vit) — [£(*)) »K*)]- Therefore, the 
variational equation of the intrinsically constrained system (M — >• Q/G, VJ 1 , / + C M ) is satisfied 
for arbitrary vertical variations in (y(i),£(t)) (i.e. vertical w.r.t the bundle M — > Q/G). 

Horizontal variations. An arbitrary horizontal variation of q(t) determines, after reduction to 
M, an arbitrary horizontal variation of y(t) in Q/G^. The corresponding variation of £(£) is not 
horizontal w.r.t the associated connection on g but the contribution of the vertical part to the 
variational equation will vanish since F/RJ* = 0. We then have £C(lZ f *) h = — ((C M )' 1 + /) which 
is 

££(W) h = i x W(y)-f(x,£) 

due to formula ([5]). 

(ii) The converse statement is easily shown in the following way. Let q(t) denote a curve in Q 
that projects onto (x(t),y(t),£(t)), the latter curve being a a critical curve of the intrinsically 
constrained system on Q /G^ x q. This implies that 



W(x, y, £)dt + / </(±, 0, Sx) + (e(y), Sy)dt = 
i Ji 

for an arbitrary variation of the critical curve. If we consider vertical variations to (y(t), £(£)) for 
which Sy(t) — 0, then F^IZ 11 = ji (i, £) — fi(y) = 0. In terms of the curve upstairs, this is precisely 

JL(q(t)) = \x. Assume that we restrict the class of variations of (y(t),t;(t)) to the projected 
variations of q(t) living in the total space Q. By using the invariance of R^, it is now easily seen 
from the previous variational equation for W 1 that 

Rf t (q)dt + f{(F + G^)(q),Sq)dt = 0. 
i Ji 

Proposition [T21 concludes the proof. 

In the previous Theorem we introduced the reduction of a Lagrangian system (Q, L, F) towards 
an intrinsically constrained Lagrangian system on a bundle M over Q/G. It should be clear that 
due to the product structure of M , there are two intrinsic constraints: F^IZ^ = or ji(x, £) = p,(y) 
and dyW 1 + C^ivY = 0. The latter constraint belongs to the gyroscopic class (see section l2~Tj) 
and in the following Proposition we study this constraint in more detail. 

Proposition 13 Let (v x ,y,£) be arbitrary in Tm{Q/G) = T(Q/G) x Q/G^ x q and (y,fj) in 
M = Q/G^ x q with f) arbitrary, then 

(d y n^(v x ,y,i),(y,fj)+Q ll ) = -{ad%fj,(y),f}). 
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Proof. The proposition follows from the following computations, with fj — [q,g(Q)] G g arbitrary 
for some curve g(e) in G through the identity at e = 0: 

(d v 1l»,(y,fj) + Q») = -{dy(ji(y),i},(3j,rj)+av) 

d 



de 



[q, (v, A d g (e)0] = (ad^A)O = -(ad|/u,J?)- 

e=0 



To provide a complete insight into the intrinsic constraint, we now compute (C)' u - Let v y G 
TylQ/Gfj,) be arbitrary, then from formula (J5J) we find 

(Oy),(2/,??)+0M) = -/ 3M K,(2/,^)+0m) = (ad^jM). 

The constraint on critical curves (y(£),£(£)) then reduces to ad^vM fl(y(t)) = ad%, t ^ jl(y(t)). Due 
to the nondegeneracy of ad*, this is equivalent to saying that y v (t) = (y(i),£(t)) + fl M . 

We have thus proven the following reduction technique. Because it preserves the Lagrangian 
nature we refer to it as Routh reduction. 

Theorem 2 Any critical curve q(t) of an invariant Lagrangian system (Q, L, F) with momentum 
fj, projects onto a critical curve (?/(£),£(£)) of the intrinsically constrained Lagrangian system 
(M -> Q/G, W, f + C M ). The Euler- Lagrange equations for (M -)■ Q/G, W , f + C) are 

££(Kn h (x,y,t) = ii&(y) ~ f(i,i) 
ji(x,i) =ji(y) 

y v = (y, £) + &»■ (9) 

4.1 Remarks on Routh reduction 

Remark 3 First we would like to draw the attention to the structure of the Routhian R^ = 
L — lu^. If we start with a critical curve of the Lagrangian system (Q,L,F) with momentum 
fj,' ^ /i (that is, a critical curve of (R^,F — G M ) with non zero momentum), the projected curve 
(y(t), £(£)) is a critical curve of a variational problem for the system (M — > Q/G, TV 1 , / + C) with 
constraints in the variations for £. The constraints are those described in Section §3.51 but w.r.t 
the action of G^. Hence, (y(£), £(£)) may not be a critical solution of the free problem defined by 
(M — > Q/G, TV 1 , f + C p )- The essential point is that, for momentum // = fi, we have ¥^71^ = 

along such projected critical curves; the constraint on St; becomes negligible and the statement 
holds for arbitrary variations of (y(t),£,(t)) in M. Moreover, the fact that ¥/JZ^ = allows us 

to separate the variations in the variables y and £. This is not possible if one works with the 
reduced lagrangian I . We could summarize this by saying that the Routhian R^ is an alternative 
Lagrangian (up to gyroscopic forces) for L such that critical curves with momentum /i become 
critical curves with zero momentum. Precisely this zero momentum condition will ensure that 
there are no constraints in the reduction of the variations. 

Remark 4 Next, we wish to explicitly mention the case where G is Abelian (for instance the 
case with cyclic coordinates). This assumption simplifies matters significantly. Firstly, we have 
that G M = G, q = Qfj_. This implies that the map 7f M is the identity and that the constraint 
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V" = (z/)O+0m m t ne Routhian reduced system (M — > QjG.W 1 , f + C 1 ) is trivial. The equations 
of motion for this system then become: 

Remark 5 In order to relate the results presented here with the existing literature, we give an 
outline of what happens when the action of G on Q is a left action. The structure of the Routhian 
reduced equations remains the same, except that some objects will change sign due to the fact 
that the structure equation now takes the form dui = fl + [u,u]. This will affect the reduced 
objects in play, for instance the 2-form /3 M becomes 

^ = (£F,ad*/I). 

The Routhian reduced equations of motion are 

ji(±,i) = jx{y) 
dyW = ad^(zj). 

Note that also d y TZ^ also changes sign in comparison with the right invariant case, so the third 
equation is equivalent to y" = (y, £) + g M . 



[T M (Q/G),T11J + ^ 



TQ/G„ 





T M (Q/G),1l%,f + C 



Figure 2: Different realizations of TQ/G^ 



Remark 6 The choice of the connection ui will alter the Routhian and force term in the re- 
duced Lagrangian system in such a way that the critical curves of the intrinsically constrained 
Lagrangian systems coincide. Assume that two principal connections ui\ and uj^ are chosen. Then 
from standard connection theory both connections are equal up to an equivariant g-valued one- 
form 6, i.e. ll>2 = Wi + 8, with R*S = Ad g -i ■ 5 and (5, c(£)) = for all £ G g. The one-form 
S reduces to a g-valued one-form on Q/G, denoted by 5. Consider the one-form 8^ = (/}, S) on 
Q/Gfj, (in fact it is a map from Q/G M to T*(Q/G)). We now study the effect on the reduction 



process. First, the connections give two possible identifications TQ/G^ 
The map </>i2 (see Fig. [2]) assumes the following form 

<f>i2(v x ,y,£) = (v x ,y,i+S(vx)). 



T{Q/G) x Q/G^ x 



If we consider the Routhian functions 1Z%, 71% and the force terms £f , £2 associated to uii and 
W2 respectively, then 



(4>* 12 n%){v xl y,C) 
tt(vy) = Cf(%) 4 



fci{v x ,y,Q 

v„dS". 



(Hy),v x ) 
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This shows that both (intrinsically constrained) Lagrangian systems are equivalent if we remem- 
ber the fact that a linearly velocity dependent potential in a Lagrangian is equivalent to the 
addition of a gyroscopic force term (see [T7]). The choice of the connection can be used to sim- 
plify the Routhian (paying the price that a gyroscopic force term is added to the system) or 
to simplify the gyroscopic force terms (paying the price that the Routhian will contain linear 
terms in the velocity). We return to this fact when we consider classical Routhian reduction 
and compare some recent geometric results in the literature and the formulation of the result as 
stated in [T7] . 

Remark 7 We end this section with a presymplectic formulation of the Routh reduced system. 
Using the notations from Section 12.11 we can consider the pull-back of the canonical symplectic 
form uiq/g from T*(Q/G) to Tm(Q/G) using the Legendre transform of W" . Since the intrinsic 
constraint associated to the fibration Q/G^ — > Q/G comes from a regular gyroscopic term, the 
following theorem is a direct consequence of the definitions in Section 12.11 Here we use the 
projections Wi : T* M {Q/G) -> T*{Q/G), ir 2 : M ->• Q/G^. 

Theorem 3 The Routhian reduced system of a Lagrangian system (Q, L, F) can be formulated 
into a presymplectic system on TQ/G^ = Tm(Q /G), where the presymplectic two-form is w-r.h = 
(tFi o WiTZ^Yuiq/q + ir^P^. The critical curves (y(t),^(t)) of the reduced system satisfy, with 

7 (t) = (i(t), !/(«),£(*)) 

(ijUJ-ilx = -dE n ^ + /)| 7 . 

In the case that (Q,L,F) is conservative (i.e. F — 0), we can apply the constraint algorithm 
to show that there exists a submanifold S of the final constraint submanifold on which a unique 
second-order vector field T exists such that 

(irUJizf = —dEiz^)\ s . 

The integral curves of T are lifted critical curves of the Routhian reduced system. 



5 Reconstruction 

The reconstruction process deals with the problem of finding a critical curve q(t) of the Lagrangian 
system (Q, L, F) such that q(t) projects onto a given critical curve (y(t), £(£)) of the intrinsically 
constrained Lagrangian system (M — > Q/G, 7^,/ + C M ) with M = Q/G^ xg-> Q/G. Using 
Theorem [TJ it is sufficient to construct a curve q(t) in Q that projects onto (y(t),£(t)) for all t, 
i.e. such that (i) [q{t)] Gft = y{t), and (ii) |(t) = [q{t),u} q{t) (q(t))] G . 

To find such a curve q(t) in Q, we choose an arbitrary starting point q a £ Q such that TT fl (q a ) = 
y{a). Next, consider the unique horizontal curve qn{t) through q a that projects onto x(t). If a 
curve q(t) exists for which (i) and (ii) holds, then it is a gauge of qh{t), that is, q(t) = qu(t)g{t) 
for a certain curve g(t) in G The reconstruction process then only consists of solving a first 
order differential equation on G to determine g(t). For that purpose, let £(£) be the curve in g 
such that [qh(t),£(t)]G = £(t) <= 3- Then, condition (ii) is satisfied if and only if g(t) solves the 
first order differential equation 

with g(a) — e. This determining condition for g{t) comes from the fact that the projection of 
q(t) onto q has to be £(£). 
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It now remains to check that (i) is true for the curve q(t) = qh{t)g{t), i.e. y(t) = [qh(t)g{t)\G„- 
Fix a time t, then the horizontal part of the tangent vector to the curve [qh(t)g(t)]G is the 
horizontal lift of x(t) to the point y(t). Similarly it follows that the vertical component of the 
curve [qh(t)g(t)]G^ is precisely (y, £) + fl M . Since g(a) = e the initial condition for [q(i)]G M is 
y(a), we may conclude that y(t) and [<?(£)]g,, coincide (both curves have the same horizontal and 
vertical parts and the same initial conditions). In this sense, the intrinsic constraint associated to 
the gyroscopic force £ M guarantees that {y{t),i;(t)) is the projection of a lifted curve q(t) in TQ, 
i.e. it reflects the second-order nature of the original Lagrangian system (Q,L,F). To conclude 
this section we unite Theorems [TJ [2] and [3] 

Theorem 4 Any critical curve of the Lagrangian system (Q,L,F) on Q with momentum fi 
projects in M = Q/G^ x g onto a critical curve of the intrinsically constrained Lagrangian 
system (M -> Q/G,H' i ,f + C M )- Conversely, any critical curve of (M -)■ Q/G,K>*,f + C p ) is 
the projection of a critical curve of (Q,L,F). The Euler- Lagrange equations of motion for the 
intrinsically constrained system are 

££(TZn h (x,y,t) = i*&(y) - /(z,f) 

y v = (y,t) + &»- 

These equations are equivalently formulated in a presymplectic setting. Let^/(t) = (x(t), y(t),£(t)) 
in Tm(Q/G) be associated to a curve (y(t),£(t)) in M. This curve is critical if and only if "f(t) 
solves 

(i^uj-n* = -dE n „ + /)| 7 . 

6 Routhian reduced systems: standard cases and examples 

6.1 The regular case 

The classical geometric description of Routh reduction (for example [T], [15]) mainly differs from 
the one described above in that the unreduced Routhian is only defined in j£ (/i) thus providing 
a reduced Routhian independent of the variable £. This is done for T — V Lagrangians. This 
classical setting is a special case of Theorem H] when additional regularity conditions are assumed 
(the constraint F1Z 11 — is a regular configuration constraint). The idea is, using the results 
from i J2.2[ to eliminate the variable £ by means of the momentum equation ji(v x , £) = p.(y). For 
that purpose we now assume that the Lagrangian L is both G-invariant and G-regular, i.e., the 
mapping Fd is invertible (see ij3. 21 above). Let m : T(Q/G) xgM T(Q/G) x g be the inverse 
ofF|L 

There is a subtlety in the application of the result from £12.21 Until now we have applied the results 
from iJ2]to Routh reduction by identifying M with Q/G^ x q and N with Q/G. In the following 
we deviate from this and we let N be Q/G^. In other words, we reinterpret the constrained 
system (M -> Q/G, 7^, / + C) as a constrained system (M ->■ Q/G^ll^, f + C M ), i.e. for now 
we forget the constraints in y defined from the fact that TV 1 depends on y but not on y. Then, 
from £12.21 it follows that the intrinsic constraint Fe7?/ i = is a regular configuration constraint, 
and using the inverse ki we find a section 7 : T{Q/G ll ) — > Q/G^ x § in the sense of £|2. 21 given by 

j(v y ) = (y,Kf(v x ,p,(y))), (10) 
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where v x — Tn^^Vy) with W^ : Q/G^ — > Q/G the natural projection; and nf — pr 2 o Ki with 
pr 2 : T(Q/G)xg — > g the projection onto the second factor. As in £|2.2[ we are now able to consider 
the Lagrangian system on N — Q/G^. However, recall that we temporarily forgot the intrinsic 
constraints in the variable y. We can reintroduce these constraints by making the remark that the 
resulting system on N carries intrinsically constraints on the hbration N = Q/G^ —> Q/G. For 
notational convenience we directly introduce the Lagrangian system on N as being a constrained 
system. This systems reads {Q/G^ — > Q/G, 1Z , / + C M ), where 

T?- : T Q/G ^Q/G) ^M : (v x ,y) *■ n^(v x ,y, K f (v x ,fi(y))) 
1 ■ T Q/G ^Q/G) -> T*(Q/G) : (v X) y) ^ /(«,(««, £(»))) 

and £** as usual the gyroscopic force associated to the 2-form /3 M on Q/G^. We conclude that the 
critical curves of the constrained Lagrangian system (M — > Q/G, W 1 , / + C M ) ai 'e in a one-to-one 
correspondence to the critical curves of {Q/G^ — > Q/G, 1Z , / + C M )- 

Recall that /3 M can be decomposed into a horizontal and vertical part: (fi, —ad*). Similar to 
previous definitions, we let d y 1Z be the shorthand notation for the restriction of the Euler- 
Lagrange operator ££(7Z ) to vectors vertical to Q/G^ — > Q/G. From Theorem [4] we easily 
have: 

Theorem 5 Any critical curve of the Lagrangian system (Q, L, F) with momentum p projects 
in Q/G^ onto a critical curve of the intrinsically constrained Lagrangian system (Q/G^ — > 
Q/G,1Z , f + C M ), and vice versa. The Euler- Lagrange eguations of motion for the intrinsically 
constrained system are 

E£(jl ,i ) h (x,y) = i x &(y)-J(x), 

d y 'Jl ll {x,y) = -ad*„fj(y) 

The equation d y TZ (x, y) = —ad yV p(y) can be rewritten as 

y v = (j/,K?(i,A(j/)))+flM- ( U ) 

This theorem is similar to the classical geometric formulation of Routh reduction as in e.g. [51115) 
where the variational problem is confined to the set J£ (p). It is known that critical curves 
in TQ with conserved momentum p are in one-to-one correspondence to critical curves of the 
variational problem under the constraint J£ (p) C TQ. In these references the Routhian is 
defined as the restriction of R^ 1 to J i ^ 1 (/i), and since this restriction R >1 '\j-^i, i \ is invariant with 

respect to the action of G M , it is reducible to J^ 1 (p)/G ll . Using the fact that in the G-regular 
case (cf. Proposition [5]) 

j£ 1 (p;)/G IM ~T(Q/G)xQ/G li , 

it is straightforward to show that this reduced Routhian is precisely 1Z . 

In 15 (as well as in the original work of Routh) the Lagrangian L is of the type L = T 2 — V with 
T 2 = 7}p(v q ,v q ) a kinetic energy associated with a non-degenerate positively definite invariant 
metric p on Q and V : Q — > R an invariant potential energy. Following [2l [9] one chooses the 
principal connection to be the mechanical connection on Q (the connection whose horizontal 
distribution is orthogonal to Vtt w.r.t the kinetic energy metric p). The induced metric on Q/G 
is denoted by p (defined by p(v x ,v x ) = p(v x ,v x )) and the vertical part of the metric defines a 
metric / on g. This is called the inertia tensor and we assume it to be regular (i.e. the Lagrangian 
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is G-regular). The momentum map then equals ji(x,£) = bj(a;)(£) = fl(y) and the Routhian on 
Q/G M is 

T^{v x ,y) = -p(v x ,v x ) - (V(x) + -(A(j/)»ttf(*)(M(y)))J • 

One should note that this Routhian may differ from the one found in classical books dealing with 
Routhian reduction. The difference lies in the fact that one typically works in a local coordinate 
chart, where the choice of the connection to is taken to be the local connection with vanishing 
connection coefficients (and thus zero curvature). In the reduced system the curvature force- 
term vanishes. This choice of connection with vanishing curvature will imply that the Routhian 
contains linear terms in the velocity, see also Remark H>] 

For the reconstruction process in this regular situation one may proceed as in fJS] Indeed, given 
a solution y(t) of the reduced problem defined by 1Z we consider the curve j{y{t)) in Q/G^ x §, 
where 7 is as in (1101) . Then a family of critical curves q(t) of the problem defined by L can be 
recovered as in the general case. Alternatively (see [HI §IV]), the curve q(t) can be recovered by 
considering any curve q(t) with y(t) = [<?(£)] g„ and writing q(t) — q(t)g(t), where g(t) £ G M is 
obtained by requiring that q(t) £ J£ (/u),Vt. 

6.1.1 The (pre)-symplectic nature of classical Routhian reduction 

It was formulated in [15] how the symplectic bundle picture for Routhian reduction looks like. 
Using Proposition Q] we can extend this to the more general case of presymplectic geometry. 
The Routhian reduced equations of motion in the classical case (so with a positive definite 
kinetic energy) are intrinsically constrained Euler-Lagrange equations living on Tqiq l (Q/G) = 
T(Q/G) x Q/Gp. Using the notations from SectionHwe have M = Q/G^, N = Q/G. We can 
say that critical curves of the Routh reduced system are solutions to the following presymplectic 
system, with 7(4) = (x(t),y(tj): 

(iKt)((7ri oFi7T)W +7^0") = -dE^, +/) | 7 , 

where F^ : T(Q/G) x Q/G M -> T*(Q/G) x Q/G^ stands for the fiber derivative along the 
T(Q/G) factor only. In the case that FiTZ is a diffeomorphism then the presymplectic form is 
symplectic. This is the case if L is a regular Lagrangian, see also [12] . 

6.1.2 The modified Tippe Top 

The purpose of this example is to demonstrate the relevance of Routh reduction for noncon- 
servative systems. In [3 a Tippe Top was described as a Lagrangian system with a dissipative 
force term. The system is invariant under the action of S* 1 and a local stability analysis for 
relative equilibria was proven by means of Routh reduction. We only mention here the setup of 
this particular Lagrangian system within theoretical framework developed in this article. The 
configuration space is Q = 5*0(3) and we will work in local coordinate chart determined by the 
Euler angles (cp, 6, ip). We consider the Abelian group G = S 1 and the action on 5*0(3) is defined 
by, for a £ S 1 : 

(ip, 6, V) ^ (<p + Rot, 0,ip — ecc), 

for e, R > 0. Next, we define the Lagrangian L which is of mechanical type and a dissipative 
force term F which represents the friction due to the slipping of the contact point of the Tippe 
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Top with the horizontal plane on which it evolves: 

L = - ((e 2 m sin 2 9 + A)9 2 + A sin 2 9>p 2 + C{ip + ip cos 9) 2 ) - mg(R - e cos 9), 
F = -fi(R- e cos 9) 2 9d9- fie sin 2 9 (op + Rip) dip - fiR sin 2 9 (Rip + vp) dip, 

It is straightforward to see that (Q, L, F) is invariant. The associated momentum map: TQ —> R 

is: 

J L = FL(q,q)(Rd<f,-ed^) = RA sin 2 9ip + C{ip + cos 9p)(R cos 9 - e) 
= R ( {A sin 2 9 + C cos 2 9) <p -Ceif) + C cos 9(Rip - tip), 

also called the Jellet integral. Fix a value \x € R of Jj,. The group is abelian, therefore G M = S 1 , 
and SO^/^ 1 = S" 2 . In order to define the Routhian, we consider the following splitting of 
TSO(3): the horizontal subspace is defined by the (invariant) subspace spanned by 

dg and ed v + Rd^. 

The curvature is zero, and we know that the Routh reduced system is a Lagrangian system on S 2 , 
i.e. of the form ££(1Z ) = — /, cf. Theorem [5] with f2 = 0. In this case, the Routhian becomes: 

R^ = L-^-^(R^-eiP) 

The next step is to use the Jellet integral and to compute R . Due to the complexity of the 
computations, we refer the reader to |3]. 

6.1.3 The free rigid body 

We discuss the example of a free rigid body to illustrate the meaning of the equation d y TZ = 
— ad*„/i(j/) in Theorem [5j We have shown that it corresponds to y v = (y,£) + g^, with £ 
determined from ji(i,£) — fj-(y) and here we illustrate this correspondence explicitly. 

Following [15] , the configuration space of the free rigid body is the entire group Q = G = SO(3) on 
which G acts from the left (in view of remark[S]we should change the sign of the 2-form: +ad*). An 
element in SO(3) is denoted by A and corresponds to the rotation taking a reference configuration 
(with principle inertia axis) of the rigid body to its configuration at time t. The Lagrangian is 
the left invariant kinetic energy and the momentum equation corresponds to the spatial angular 
momentum. Using the notations from |15j . we assume that the the fixed momentum equation 
corresponds to 7r = /xk (where w is the spatial angular momentum and k the unit vector in 
the positive z-axis). We choose the standard connection ui(A)(A) = AA^ 1 , and we identify the 
(trivial) bundle g with g by [A, £]g h- A^ 1 S,A, i.e. the representant of L4,£]g at unity. The 
map determined by the connection u which is used to identify TSO(3)/G with T(Q/G) x g = g 
becomes (A, A) H> A~ X A, i.e. the image is the angular velocity in the body reference frame. In 
a local coordinate neighborhood determined by Euler angles, the projection TSO(3) — > g = M 3 
equals 

with 

£.1 = 9 cos ip + (p sin 9 sin ip 

£2 = — 9 sin ip + (p sin 9 cos ip 

£3 = ip + <pcos9. (12) 
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The reduced Lagrangian I then determines a function on g and equals I — \{IiS, 2 + h£,2 + ^3^1)) 
with (Ii, 12,1s) the inertia tensor in the body reference frame. 

The isotropy subgroup G M = S 1 consists of the rotations about the k-axis; therefore ir^ : Q — > 
Q/Gfj, — S 2 : A n- A~ l \a. In terms of the Euler angles, this projection is (<fi, 9, ip) n- (8, i/j). The 
map fl : Q/G^ — > g* = g* is the momentum /ik expressed in body coordinates: 

(9,ip) | — > (fi sin 9 sin ip, fi sin 9 cos ip, [i cos 9). 

The configuration space of the Routh reduced system is simply Q/G^ — S 2 and the only mean- 
ingful equation from Theorem [5] is d y lZ = +ad* 1 ,/i(y). 

In the following we will first compute the equation £ = K 2 (ji(y)) (see (jTDJ)) and study its projection 
to the tangent bundle of Q/G^ — S 2 . Secondly, we show that this projection corresponds to the 
equation d y 1Z — +ad*i,/x(y). 

To compute k 2 we start with the reduced momentum equation j/(£) = ji{y)- 

Ii€i = /J.sin8sinip, 

hS,2 = /J,sin9cosip, 

hS,3 = MCos6». (13) 

It should be clear that since Ii ^ 0, we can compute £ explicitly 

£i = — fisin8sinip, £2 = y- fJ, sin9 cos ip, £3 = —ficos9. (14) 

Jl -<2 J3 

Next, we need to compute the corresponding element in TS 2 . In the general theory this is done 
by identifying TS 2 with the quotient 7f*g/g M . In terms of the Euler-angles, this boils down to 
computing 9 and ip out of ([15)1 : = £1 cos t/> — £2 sin t/> and ^ = £3 sin — (£1 sin ip + £ 2 cos ■0) cos 0. 
If we substitute the values for £1,2,3 in these two expressions, we find the equation y" = (y, £)+0 M , 
with £ determined from jz(£) = jj.{y)'- 

9 = fi sin sin t/> cos ip I — — 

\ * 1 -* 2 

/ 1 /sin 2 -0 cos 2 V 1 

The second way to retrieve these equations is using the Routhian 7?. . By definition it is obtained 
from W* = 1(g) - (fi(y),g), where £ is substituted by (fT4|) : 

-n /n M /x 2 / sin 2 sin 2 ^ sin 2 # cos 2 ^ 



#W) = -^- 



2 \ h h h 

Next, we compute the 2-form ad* = /3 M on TS 2 . Recall that it is the projection of dw 1 * to 
Q/Gfj,. In Euler angles, we have that (fl, u) (A, A) — jx(y) ■ (£) T = /j,(<f> + cosOip). Therefore 
/3'' = —/ism9d9 A dip, the gyroscopic force term is £'' = ~iy^ L — —fJ,sin0ipd9 + [isin99dip and 
the Euler-Lagrange equations for the Routh reduced system are ££(1Z ) = — £ M : 

d ^ 2 • 2 a ■ , 1 ( 1 l \ an 

-fi sin 9 sin ip cos ip I = —fismOtf 



dip \Ii h 

dlt 2 . n n ((sm % ip cos 2 V\ 1 1 

= — fi sin w cos w — 1 — — — I = /j, sin 9ip. 



89 " VV h h J h, 

It should be clear that these Euler-Lagrange equations, when brought into normal form, are 



precisely the equations of motion obtained above by projecting £ = K 2 (p(y)) to TS 2 . 
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6.2 Lagrangian systems subjected to magnetic forces 

Assume that we study a conservative Lagrangian system L, where L = T2+T1 — V, T2 — \p(v q , v q ) 
is the kinetic energy and T\ — (a(q),v q ) with a and invariant 1-form. These linear terms in the 
kinetic energy are found, for example, when the mechanical system is subject to magnetic force 
terms (for instance, charged particle/rigid body evolving in a magnetic field). The restriction 
of a to horizontal vectors is denoted by a h and is projectable to a one-form a on Q/G. The 
restriction of a to verticals determines a section a v of g* — > Q/G. We consider the mechanical 
connection u> defined by the metric p as before, and let / be a regular inertia-tensor (i.e. L is 
G-regular). The momentum condition then reads 

with inverse £ = tt/(/i — a"). The Routhian on Q/G^ assumes the form 

T^{v x ,y) = 2~P~{v x ,v x ) + (a(x),v x ) - (V(x) + -((p. - & v )4j(P>- a v ))j . 

6.2.1 The heavy top in a magnetic field 

We consider a constant magnetic field along the z-axis in the inertia frame. The Lagrangian of 
a charged rigid body with a fixed point is invariant under the left action of rotations about the 
z-axis (i.e. when (<j), 9, ip) is a coordinate chart determined by the Euler-angles, then <ft is (locally) 
cyclic). Thus Q = 50(3), G M = G = S 1 and Q/G = 5 2 , see also [15]. The Lagrangian equals, 
with £, as written down in Eq. (|12j) : 



1 - 



2 

4=1 

with f2 = q/mc, with q the charge of the heavy top, g the earth's acceleration, m the mass and x 
the vector with length e connecting the fixed point with the center of mass in the body reference 
frame. We assume that the top is symmetric and that the moving frame is directed along its 
symmetry axis (i.e. I\ — h). In this case, the Lagrangian becomes 

L = -hd 2 + -h"^ 2 +I 3 cos9<j)tp + -p^ 2 - mge cos 9 - QBip^ + I 3 cos 8ip), 

with P00 = I\ sin 2 9 + I 3 cos 2 9. Thus a = —^lB{p t j }t j } d<p + ^3 cos9dip), a v = —ClBp^ and a h = 0. 
Here we have chosen uj to be the mechanical connection u = d<j)+{Is cos 9dip)/ p^. The conserved 
momentum Jl reads 

^3 cos 9ip — flBptjttj, = p. 



The next step is the computation of the Routhian 1Z , which can be carried out in two ways: 
either one computes <fr out of the momentum equation, and then substitute the expression in 
L — p((f> + (^3 cosOip)/ ptfxf/), or one computes p and use the expression of 1Z determined above. 
After straightforward but tedious computations one finds that the Routhian 1Z for the magnetic 
heavy top assumes the following form 
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with 

2" 



V li (O)=mgcos0+-- 



In the case that FIB C 1, we may neglect the term in (J7B) 2 and then the effective potential 
Vfj, equals up to constant the effective potential for the heavy top (i.e. set B = 0). In this case 
the critical curves of both problems coincide, but the reconstruction process will differ. The 
reconstruction equation, roughly said, is the inverse of the momentum map: 

P4>4> 

This result is also known as Larmor's theorem (for example, see [5]): the motion of a charged 
particle in a constant magnetic field is the motion of the particle in the absence of the magnetic 
field superposed with a rotation with constant frequency QB, the Larmor frequency. 

6.3 Lagrangians linear in the symmetry generators 

Here we study the case where the Lagrangian of an invariant conservative Lagrangian system 
(Q, L, 0) generates a reduced momentum map ji that does not depend on the variable £, i.e. for 
a critical curve (y(t),£(t)) projecting onto x(t): 

jl(x) = p(y). 

This situation occurs for instance if one studies Lagrangian systems for which the kinetic energy 
is not strictly positive definite, i.e. if L(v q ) = ip(q)(v q , v q ) — V(q), with p and V right invariant, 
p a symmetric two-tensor such that p(cr(£), c(£)) = for any £ S Q. 

In this situation we can consider the metric p on Q/G defined in the usual way as ~p(v x ,v x ) = 
p(q)((v x )q , (v x ) q ), for ir(q) — x and v x arbitrary in T X (Q/G). The remaining part of p yields a 
0*-valued one-form p on the quotient space Q/G. The momentum constraint for a critical curve 
then yields p(x)(x) — p{y). 

We can apply Theorem [5] and we find that the Routhian 1Z^(v x ,y,£) equals 

-p(x)(v x ,v x ) - V(x) + ((p(x) - jS(y)),|). 

From Section [2l Proposition [2] we know that the critical curves of this intrinsically constrained 
Lagrangian system (M — » Q/GjTZ^X^) project onto critical curves of the Lagrangian system 
(Q/G p, 1Z , C M ) constrained to the submanifold in T(Q/G M ) determined by p(x)(v x ) = p,(y). The 
function 1Z is defined as 

T^(v x ) = -p(x)(v x ,v x ) - V(x), 

where we identified V with its reduced function on Q/G. This particular situation is encountered 
in general relativity, where solutions to Einstein's equations admit lightlikc Killing vectors. 

6.3.1 Lightlike Killing vectors 

This example shows how the presented theory can be applied to reduce geodesic equations in 
general relativity where the metric admits a lightlike Killing vector. We assume that the Killing 
vector X is complete, i.e. the flow of X determines an Abelian group action of R on the spacetime 
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(M, g). Since X is Killing it determines a symmetry of the kinetic energy Lagrangian L(v q ) — 
\g{v q ,v q ). 

A relevant local example from general relativity illustrating the above can be found in the so 
called pp-wave solution to Einstein's equations (see [TT]): 

g = H(u, x, y)du 2 + 2dudv + dx 2 + dy 2 , 

expressed in the Brinckmann coordinates. The vector field d v is Killing, and the momentum 
constraint then assumes the form Jl(u, x, y,u) = u = (i. The geodesies with momentum /x are 
critical curves for the variational system on TR 3 with Lagrangian equal to 

n tl = ^( f i 2 H(u,x,y) + x 2 + y 2 ) 
constrained to the submanifold u = \i. 
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